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Abstract 

We apply a quadratic hedging scheme developed by FoUmer, Schweizer, and Sondermann to 
European contingent products whose underlying asset is modeled using a GARCH process and show 
that local risk-minimizing strategies with respect to the physical measure do exist, even though an 
associated minimal martingale measure is only available in the presence of bounded innovations. 
More importantly, since those local risk-minimizing strategies are in general convoluted and difficult 
to evaluate, we introduce Girsanov-like risk-neutral measures for the log-prices that yield more 
tractable and useful results. Regarding this subject, we focus on GARCH time series models with 
Gaussian innovations and we provide specific sufficient conditions that have to do with the finiteness 
of the kurtosis, under which those martingale measures are appropriate in the context of quadratic 
hedging. When this equivalent martingale measure is adapted to the price representation we are 
able to recover out of it the classical pricing formulas of Duan D95' and Heston-Nandi [HNOO] . as 
well as hedging schemes that improve the performance of those proposed in the literature. 

1 Introduction 

GARCH models [E82[ 1686) IPGE93' have been introduced in the modehng of the time series obtained 
from financial stock prices with the objective of capturing via a parametric and parsimonious family of 
processes several features that have been empirically documented and that escape to more elementary 
modeling tools. For example, the constant variance and drift time series model that one obtains out of 
the strong Euler discretization of the lognormal model that underlies the Black, Merton, Scholes (BMS) 
option valuation formula |BS72[ IM76| is not able to account neither for the volatility clustering in the 
time series of the associated returns nor for the leptokurtosis (fat tails) in their distribution. Moreover, 
the oversimplification in modeling the stock returns is a source for the appearance of contradictions in 
the implications of the BMS pricing formula, like the smile shaped curve that one observes when the 
implied or implicit volatility is plotted as a function of either moneyness or maturity. 

From the modeling point of view, the GARCH family is successful at the time of reproducing the 
above mentioned empirically observed features. Moreover, these models are particularly attractive 
from the mathematical point of view since the conditions for the existence of stationary solutions can 
be simply formulated and, additionally, most of the standard techniques in the time series literature 
concerning model selection and calibration can be adapted to them (see for instance |G97[ IH94| and 
other standard references therein). 

The situation becomes more complicated when we try to price contingent products whose underlying 
asset is assumed to be a realization of a GARCH process. The discrete time character of the model, 
together with the infinite states space usually assumed on the innovations, makes the associated market 
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automatically incomplete, in the sense that there are payoffs that cannot be replicated using a self- 
financing portfolio made only out of a bond and the risky asset. This difhculty has been extensively 
treated in the literature using various approaches. 

A first way to address this problem (see |D95l IHNOOj ) consists of adding a term to the GARCH 
model in the spirit of the NGARCH and VGARCH models introduced in |EN93| ; when the conditional 
mean of this modified model is computed using the physical probability, a term proportional to the 
conditional variance appears that makes it non-risk-neutral. The associated constant of proportionality 
is interpreted as a return premium per unit of risk. The price of the options that have this model as 
underlying is then defined as the discounted expected value of the payoff with respect to a new pricing 
measure under which the GARCH process is risk-neutral; other requirements are also imposed on the 
pricing measure to ensure, for instance, that the Gaussian character of the innovations is preserved 
under the new measure. In the case of jD95| . a utility maximization argument gives legitimacy to this 
definition. 

A different approach consists of finding continuous time processes that extend GARCH to that 
setup, following a scheme that makes the market complete and yields a price formula in the spirit of the 
Black-Merton-Scholes theory (see [KT981 IKV081 |D97] and references therein). It is worth mentioning 
that tackling the problem in this way, Kallsen and Taqqu |KT98j obtain results that are consistent 
with those of |D95| as far as the pricing formulas is concerned, but disagree on the associated hedging 
strategies (see |GR98| for a discussion). 

In this paper we will focus on the hedging side of the problem and will implement in the GARCH 
setting the quadratic hedging approach developed by FoUmer, Schweizer, and Sondermann (see |FS86[ 
IFSc901 ISchQlj , and references therein) . Given a probability measure, the theory developed in the papers 
that we just quoted gives a prescription on the construction of a generalized trading strategy that 
minimizes the local quadratic hedging risk (to be defined later on). Quadratic hedging techniques 
can be subjected to improvement since they do not make a difference between hedging shortfalls and 
windfalls, which should be obviously treated differently as far as the associated risks is concerned. 
Even though this point has been addressed in a variety of works (see [POOj . and references therein) 
the associated hedging and pricing problem is more convoluted; we will hence put off the use of these 
techniques in the GARCH context to a future work. 

The contents of the paper are organized as follows: Section [2] contains a quick review of the GARCH 
models as well as the notions on quadratic hedging that are used later on in the paper. The last part 
of this section contains a first result that shows the availability of the quadratic hedging scheme in 
the GARCH context and a second one where we spell out the conditions under which there exists a 
minimal martingale measure; whenever this measure exists, the value process of the local risk-minimizing 
strategy (with respect to the physical measure) admits an interpretation as an arbitrage-free price for 
the derivative product we are dealing with. Unfortunately, the range of situations in which the minimal 
martingale measure exists is rather limited and, as we will see, is constrained to GARCH models 
with bounded innovations; this limitation is, from the modeling point of view, not always appropriate. 
Moreover, the expressions that determine the optimal hedging strategy using this measure are convoluted 
and hence of limited practical applicability. 

The situation that we just described motivates us to carry out in Section [3] the local risk mini- 
mization program for a well chosen Girsanov-like equivalent martingale measure. We implement this 
program for GARCH models whose innovations are Gaussian. Quadratic hedging with respect to a 
martingale measure yields much simpler expressions, admits a clear arbitrage-free pricing interpretation 
and, additionally, the corresponding strategies minimize not only the local risk and the quadratic risk, 
but also the so-called remaining conditional risk (these concepts will be defined later on). Moreover, we 
will prove that a linear Taylor expansion in the drift term of the local risk minimizing value process with 
respect to this martingale measure coincides with the same expansion calculated with respect to the 
physical measure; consequently, since in most cases the drift term is very small, one can safely compute 
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the risk minimizing strategy with respect to the martingale measure, which is much more convenient, 
and one obtains practically the same value had one used the much more convoluted expressions in terms 
of the physical measure. 

Even though the equivalent martingale measure always exists, the quadratic hedging scheme requires 
the process modeling the underlying asset to be square summable with respect to the pricing measure. 



In Theorem 3.1 we show that a sufficient (but not necessary) condition for that to hold is the finiteness 
of the kurtosis with respect to the physical measure. It is worth mentioning that with this change of 
measure, the independent Gaussian innovations of the original GARCH process remain automatically 
independent and Gaussian after risk neutralization and there is no need to impose this feature as an 
additional condition (compare with, for example. Assumption 2 in |HNOO] ). 

The developments that we just explained are formulated using log-prices as the risky asset. In 



Section 3.1 we reformulate these results in the price representation using predictible drift terms. This 



degree of generality allows us to recover in Section 3.2 the classical pricing formulas of Duan |D95) 
and Heston-Nandi (HNOOj in the context of the local risk minimization scheme. It is worth mentioning 
that apart from providing an alternative interpretation to existing pricing formulas and improving some 
theoretical issues (like, for example, dropping assumptions invoked in [HNOOj regarding the Gaussian 
nature of the innovations after risk neutralization), the scheme that we propose comes together with a 
hedging scheme that performs better than the one proposed in |D95| and is simply not available in the 
case of |HNOO| . More explicitly, local risk minimizing hedging strategies with respect to a martingale 
measure minimize by construction the mean square hedging error, when this error is measured with that 
martingale measure and are hence, from this point of view, preferable to the hedging ratios proposed 
in |D95| or [KT98| . Even though there are no theoretical results that guarantee that that this difference 
in hedging performance still exists when we move to the physical probability, we have carried out 
numerical tests in Section IXSl that seem to indicate that this is indeed the case. 

Conventions and notations: all along this paper we will use the riskless asset as numeraire in 
order not to carry the riskless interest rate in our expressions. Given a filtered probability space 
(J7, P, {Tn}n£fi) and X, Y two random variables, we will denote by iJ^l-'^^] := E[X\J^n] the conditional 
expectation, cov„(X, F) := cov{X, Y\J^„) ■= En[XY] — En[X]En[Y] the conditional covariance, and by 
var„(X) := En[X'^] — _E„[X]^ the conditional variance. A discrete-time stochastic process {Xn}neri is 
predictable when Xn is J>i_i-measurable, for any n € N. 

Acknowledgments: I thank Josef Teichmann for generously spending his time in going carefully 
through the paper, giving it a serious thought, and for providing me with relevant feedback that has 
much improved it. I have also profited from discussions with Stephane Chretien, who has shared with 
me his vast culture in statistics and probability theory. Interesting discussions with Alberto Elices, 
Friedrich Hubalek, and Andreu Lazaro are also acknowledged, as well as the comments and criticisms 
of two anonymous referees that have significantly improved the paper. 



2 The GARCH family and pricing by local risk minimization 

In this section we introduce the general family of times series that we will use for the modeling of the 
stock prices. We then briefly review the basics of quadratic hedging, and we finally prove the existence 
of that kind of strategies in the GARCH context. 

2.1 The GARCH models 

Let (ri, J^, P) be a probability space and {e„}„gN ^ IID(0, 1) a sequence of zero-mean, square integrable, 
independent, and identically distributed random variables. We will denote by {J-'n}neN the filtration 
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of T generated by the elements of this family, that is, ■= cr(ei, . . . ,e„), n > 1, is the a-algebra 
generated by {ei, . . . , e„}. We will assume that J'o is made out of 51 and all the negligible events in J^. 

GARCH models were introduced in [B86| as a parsimonious generalization of the ARCH models 
used by Engle |E82| in the modeling of the dynamics of the inflation in the UK. This parametric family 
has been modified in various forms to make it suitable for the modeling of stock prices and contains an 
important number of different models. Even though the treatment that we will carry out in the rest 
of the paper is valid for all the models in the literature, we will now pick one of them, namely the one 
introduced in |DGE93| , to illustrate the features of these models that will be of relevance in the rest of 
the paper. 

An example: the asymmetric GARCH model: let {^nlngN be the sequence that describes the 
price of the risky asset that we are interested in. The asymmetric GARCH(p,q) model |DGE93[ 
IHT99) determines the dynamics of the prices {5'„}„gN by prescribing the following dynamics for the 
log-returns r„ := log(S'„/5'„_i), that amounts to a recursive relation for the log-prices s„ := log(S'„): 



P 1 



(2.1) 
(2.2) 



i=l 



i=l 



where f„ = r„ — E[rn] — r„ — fi, {e„}„gN ^ 11DN(0, 1), and u > 0, ai,l3i > 0, I7I < 1 are constant 
real coefficients. Notice that the fact of working with log-prices implies that the price process {5'„}„gN 
determined by (2.1 1 and (2.2) is always positive. The parameter 7 controls the asymmetric influence 
of shocks: if they are positive, negative past shocks raise more the variance than comparable positive 
shocks. This is an empirically observed feature of stock markets. The following proposition, whose 
proof is sketched in the appendix, characterizes the constraints on the model parameters that ensure 



the existence and uniqueness of a weakly stationary solution for (2.1|-(2.2). 



Proposition 2.1 Suppose that w > 0, at, Pi > and \j\ < 1. Then the model (2.1]-(2.2] admits a 
unique weakly (second order) stationary solution if and only if 



(H- 7^) (ai -I- ••• + ap) + /3i + •••+/?,< 1, 



in which 



var(r„) = E[a„] = ^-^^—^ 



(2.3) 
(2.4) 



Finite kurtosis: apart from the second order stationarity studied in the previous statement, it will 
be important in the developments of the paper to be able to characterize the situations in which the 
solutions of (2.1|-(2.2) have finite kurtosis, that is, they have a fourth order moment. This point is 
relevant for two reasons: first, in the presence of finite kurtosis, it can be shown that the square of a 
GARCH process is a linear ARMA model; given that for ARMA there exists a large array of preliminary 
estimation tools for model selection and calibration, one can take advantage of this situation in the 
calibration of the GARCH process that one is interested in. Second, as we will see in Theorem |3.H 
having finite kurtosis is a sufficient condition so that pricing by local risk-minimization is available with 
respect to a risk-neutral measure. In the particular case of the asymmetric GARCH model, the existence 
of finite order moments has been characterized in |LMc02bl ILLMc02] in the following compact form: 
the necessary and sufficient condition for the existence of the moment of order 2m is that 



p[E < 1, 



(2.5) 
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where (8) denotes the Kronecker product of matrices, p{B) = max {| eigenvalues of the matrix B\}, A is 
the matrix given by 



A 



OLlZt ■■■ OLpZt 


PiZt 




-^(p-l)x(p-l) 0(p_l)xl 


0(p_l)xg 




Qfl • • • ttp 






0(,-l)xp 


-^(<?-l)x(g-l) 


0(g-l)xl 



^pxp is the p X p identity matrix, and Zt := (|et| — 7et) • For m = 1, the condition (2.5 1 is the 
same as (|2.3|). The kurtosis is finite whenever (2.5) holds with m = 2. For example, in the case of a 



GARCH(1,1) model, (2.5) amounts to the following inequality relation among the model parameters: 

^2 + 2^a(l + 7^) + [(1 + ^2)2 _^ 4^2j ^ -^^ 

The paper |LMc0 2' contains the corresponding characterization for the finiteness of the kurtosis of other 
asymmetric GARCH(1,1) processes (like GJR-GARCH) or driven by non-normal innovations. 

Volatility clustering and leptokurtosis: GARCH is successful in capturing these two features that 
one empirically observes in stock market log-returns. Actually, in the GARCH context, these two notions 
are intimately related in the sense that one can say that heteroscedasticity (volatility clustering) causes 
leptokurtosis (heavy tails) and vice versa. Indeed, since we are using Gaussian innovations, we have 

This allows us to write down the kurtosis (standardized fourth moment) as 

V 2 



/C = 



5E [alel] 



3E 



{En-iHel]y 
{EKel] 



E 



3 + 3- 



{E,,^,[alel]f ~E[E,,^,[alel]\ 



{EWiei]y 



3 + 3 



5E [alel] 
var 



{E[alel]Y 



3 + 3- 



var 



{E[<yl]Y 



where E [tr^j is determined by (2.4). Notice that this expression, due to jG97j . proves that the excess 
kurtosis is positive whenever the variance of the volatility is non-zero. 

General GARCH models: the results that we will prove in this paper apply beyond time series 
models that follow exactly the functional prescription determined by expressions (2.1 ) and (2.2 ). In our 
discussion it will be enough to assume that the log-prices evolve according to: 



log 



Sn 



where {e„}„gN ^ 
to J"„_i := cr(ei 



(2.6) 

~ '''n ('''" — 1 ' ■ • ■ ' '''n— max(p,g) J ^n— 1 5 ■ • ■ : ^n— 9)7 (2-7) 

{A'nIneN is a predictable process (that is Pn is measurable with respect 



110(0,0-2 

. . ,e„_i), for all n G N), and the function a^{u. 
constructed so that the following two conditions hold: 

(GARCHl) There exists a constant w > such that cr^ > w. 
(GARCH2) The process {cr„e„}„gN is weakly (autocovariance) stationary. 



; ^ n— max 



■ : ^n—q) is 



A process {sn}„gN determined by (2.7) and (2.7) will be generically called a GARCH(p,q) process. 
Notice that (2.7) implies that the time series {cr„}raeN is predictable; this feature is the main difference 
between GARCH and the so-called stochastic volatility models. 
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2.2 Local risk minimizing strategies 

In the following paragraphs we briefly review the necessary concepts on pricing by local risk minimization 
that we will needed in the sequel. The reader is encouraged to check with Chapter 10 of the excellent 
monograph |FSc04) for a self-contained and comprehensive presentation of the subject. 

Let H(St) be a European contingent claim that depends on the terminal value of the risky asset 
Sn- In the context of an incomplete market, it will be in general impossible to replicate the payoff H by 
using a self-financing portfolio. Therefore, we introduce the notion of generalized trading strategy, 
in which the possibility of additional investment in the numeraire asset throughout the trading periods 
up to expiry time T is allowed. All the following statements are made with respect to a fixed filtered 
probability space (f2,P, J", {-^n}„e{o,...,T})- 

Definition 2.2 A generalized trading strategy is a pair of stochastic processes {£,'^,0 such that 
{£n}ne{o,...,T} adapted and {£n}ne{i....,T} predictable. The value process V of {£^,£) is defined 
as 

Vo := ^0, and Vn + ' Sn, n > 1. 

The gains process G of the generalized trading strategy (^'^,^) is given by 

n 

Go and Gn ■■=^£k ■ {Sk - Sk-i), n = 0,...,T, 

k=l 

and the cost process C is defined by the difference 

Gn '■=Vn~ Gm n = 0, . . . ,T. 

It is easy to check that the strategy is self- financing if and only if the value process takes the 

form 

n 

Fo=e? + ei-5o and Vr, = Vo + Y,^k-{Sk-Sk-i) = Vo + Gn, n = l,...,T, (2.8) 

fc=i 

or, equivalently, if 

Vo = Co = Ci = . . . = Cr. (2.9) 

Definition 2.3 Assume that both H and the {Sn}ne{n....,T} are L'^{Q,,V). A generalized trading strategy 
is called admissible for H whenever it is in L^(r2,P) and its associated value process is such that 

Vt = H, F a.s. and Vt G L^(f7,P), for each t, 

and its gain process Gt G £^(51, P), for each t. 

The next definition introduces the strategies we are interested in. 

Definition 2.4 The local risk process of an admissible strategy (^",0 process 

Rt{e,0 -.^EtiiGt+.-Gt)'], t = 0,...,T-l. 

The admissible strategy (^,^) is called local risk-minimizing if 

for all t and each admissible strategy (C'^i^)- 
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It can be shown that |FSc041 Theorem 10.9] an admissible strategy is local risk-minimizing if and 
only if the cost process is a P-martingale and it is strongly orthogonal to S, in the sense that cov„(S'„_|_i — 
Sn,Cn+i — Cn) — 0, P-a.s., for any t — 0, ... ,T — 1. An admissible strategy whose cost process is a 



P-martingale is usually referred to as mean self-financing (recall (2.9) for the reason behind this 
terminology). Once a probability measure P has been fixed, if there exists a local risk- minimizing 
strategy (^",0 with respect to it, then it is unique (see |FSc04| Proposition 10.9]) and the payoff H can 
be decomposed as (see |FSc041 Corollary 10.14]) 



H = Vo + Gj 



(2.10) 



with Gn the gains process associated to (C°,C) and L„ := C„ — Co, n = 0, ...,r. Since (C°,C) is 
local risk-minimizing, the sequence {Ln}n£{o,...,T}^ that we will call global (hedging) risk process^ 
is a square integrable martingale that is strongly orthogonal to S and that satisfies Lq = 0. The 



decomposition (2.101 and (2.8) show that Lx measures how far H is from the terminal value of the 



self-financing portfolio uniquely determined by the initial investment Vq and the trading strategy ^ 
(see |LL081 Proposition 1.1.3]). 

2.3 Local risk minimization in the GARCH context and minimal martingale 
measures 



As we pointed out in the previous section, the local risk-minimization aproach to hedging demands 
picking a particular probability measure in the problem. Given a contingent product on a GARCH driven 
risky asset, the physical probability measure is the most conspicuous one since, from the econometrics 
point of view, it is the measure naturally used to calibrate the model. 

Our next proposition shows that a local risk-minimizing strategy with respect to the physical measure 
does exist in the GARCH context. Given the specific form of (2.6) and (2.7), it is more convenient to 



reformulate the problem by finding a local risk-minimizing strategy in which we take the log-prices s„ 
as the risky asset and /i(st) := -ff (exp(sT)) as the payoff bmction. 



Proposition 2.5 Consider a market with a single risky asset that evolves in time according to a 
GARCH process satisfying (2.6) and (2.1), driven by innovations {e„}„gN ^ IID(0,(t^). Suppose also 
that the GARCH process has a hounded drift {/in},igN; that is, /i„ < i? for some _B S M and for all 
rt G N. Let h G i^(f2,P, J^x) he a contingent product on s ^ log(S'). Then, there exists a unique 
local risk-minimizing strategy for h with respect to the physical measure P, uniquely determined by the 
following recursive relations 



CT^CTfc 



fe-i 



h\l- 



1 ^ CT-l 



£,T — — Et-1 [hcr] , 



Vk — Ek 
Vt = h. 



1 I ^ 

/ill ^ CT 



1 5 er-i 



1 - 



1 - 



-Ck+1 



, k = l, 



,T- 1, 



,,T-1, 
(2.11) 
(2.12) 

(2.13) 
(2.14) 



The position on the riskless asset is given by ^ := Vk — £,kSk 



Remark 2.6 The condition on the boundcdness of the drift is verified for most models considered in 
the literature in the context of derivatives pricing. In some instances the drift is just a constant (see for 
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example |BEM07| V Other important situation where this hypothesis trivially holds is the model in |D95| . 
where the drift is given by /j,„ := r + Act„ — ^ct^; r and A are positive constants that account for the 
continuously compounded one time step risk-free interest rate and the unit risk premium, respectively. 
In this situation, it is clear that fi„ < r + ^X^ for any n e N. 



Proof of Proposition 



and (2.7) imply 



2.5 



We start by noticing that since (T„ is J^„_i-measurable, the relations (2.6) 



var„_i[s„ - s„_i] ^ cr^crl, 



(2.15) 
(2.16) 
(2.17) 
(2.18) 



for any n G { 1 , . . . , T} . 

The first fact that we need to check is that the GARCH context fits the framework established by 



Definition 2.3 to carry out hedging by local risk minimization. More explicitly, we have to verify that 
the log-prices s are square integrable. This is a consequence of hypothesis (GARCH2); indeed, for any 
n G {1, . . . , T}, s„ = So + J2i=i Mi + fulfil H 1- (^nen- Then, 



E[sl] = E 



\ i — 1 J i—1 i'Cj — 1 



Let i < j, then E [aiajSiej] = E[Ej^i[ai(Tjei€j\] ~ E[aiajeiEj-i[ej]] — E[aiaj€iE[tj\\ — 0. Since by 
hypothesis (GARCH2) E[afef] < oo and the drift is bounded, we have that 



E[sl] - E 



So 



^ E[<j^e^] < {so + nBf + ^ E[a'^4] < oo. 



i=l 



as required. 

Now, according to |FSc041 Proposition 10.10], the existence and uniqueness of a local risk-minimizing 
strategy is guaranteed as long as we can find a constant C such that (£"„_! [s„ — s„_i])^ < C-var„_i[s„ — 
Sn-i], P-a.s. for any n. In our case it suffices to take C = B"^ /{a'^uj). Indeed, with this choice and 



using (2.16) and (2.181, 



{En-l[Sn 



< 



var„_i[s„ 



= c, 



(2.19) 



as required. The recursions (2.11 )-(2.14) follow by rewriting expression (10.5) in [ FSc04j using the 
equalities ([2T5l)-([2l8|. ■ 



Expressions (2.11 )-( 2.14) are convoluted and difficult to evaluate. Moreover, expression (2.131 does 



not allow us to interpret Vk as an arbitrage free price for h at time k. There are two possibilities to go 
around this problem: the first one consists of dropping the physical probability and of choosing instead 
an equivalent martingale measure that has particularly good properties that make it a legitimate proxy 
for the original measure. This is the path that we will take in the next section. 

As an alternative, one may want to look for an equivalent martingale measure for which the value 
process of the local risk-minimizing strategy with respect to the physical measure can be interpreted as 
an arbitrage free price for h. This is the motivation for introducing the so-called minimal martin- 
gale measure. This measure is defined as a martingale measure P that is equivalent to the physical 
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probability P and satisfies the following two conditions: E 



< oo and every P-martingale 



M G L^(f2,P) that is strongly orthogonal to the price process s, is also a P-martingale. This measure 
satisfies an entropy minimizing property |S ch01, Proposition 3.6] and if E denotes the expectation with 
respect to P, then the value process Vk in (2.131 can be expressed as (see Theorem 10.22 in |FSc04] ) 



Vk = Ek[hl 

which obviously yields the interpretation that we are looking for. 

As we see in the next proposition, minimal martingale measures exist in the GARCH context only 
when the innovations are bounded (for example, when the innovations are multinomial) and certain 
inequalities among the model parameters are respected. 



Proposition 2.7 Using the same setup as in Proposition \2.5\ suppose that the innovations in the 
GARCH model are bounded, that is, there exists K > such that < K, for all k = 1,...,T, 
and that this bound is such that K < a^s/ujjB, with lu > the constant such that cr^ > w (see 
condition (GARCHl) ) and B gM. the upper bound for the drift. Then, there exists a unique minimal 
martingale measure P with respect to P. Conversely, if there exists a minimal martingale measure then 
the innovations in the model are necessarily bounded. 

Whenever the minimal martingale measure exists, its Radon-Nikodym derivative is given by 



n 

fc=i 



(2.20) 



Proof. We start by recalling that in the proof of Proposition 2.5 we showed in (2.19) the existence 
of a constant C such that 

{En_i[sn - s„-i])^ < C ■ var„_i[s„ - S„_i] 



for all t = 1, . . . , T. In view of this and Theorem 10.30 in fFSc04] , the existence and uniqueness of a 
minimal martingale measure P is guaranteed provided that the following inequality holds 



(s„ - Sn^i)En^i[Sn - Sn-l] < £^n-l[(s« " Sn-l)^]- 



(2.21) 



By (2.16) and (2.17), this inequality is equivalent to e„ < u'^Un/^-n and it obviously holds if the 
innovations are bounded and the bound satisfies K < a'^^/uj/B. Conversely, suppose that there exists 
a minimal martingale measure; Corollary 10.29 in |FSc04j implies that (2.21) holds and hence so does 
e„ < a'^Unl Hn- Given that cr„ and /i„ are J-n-i measurable and e„ is J-n measurable, this equality can 
only possibly hold whenever the innovations e„ are bounded. 



As to expression (2.20), it is a consequence of Corollary 10.29 and Theorem 10.30 in |FSc04] . Ac- 
cording to those two results, the density dP/dP is the evaluation at T of the P-martingale 



Zt (1 + Afc • [yk - Vk-i)) 



fe=i 



where Xk '■= —Ek-i[sk — Sk-i\/ysxk-i[sk — Sk~i\ and yk is the martingale part in the Doob decomposition 
of Sfc with respect to P. Therefore, we have 



Vk - Vk-i ^ Sk- Sk-i - Ek-i[sk - Sfc_i] 



Using (2.16) and (2.18) in these expressions the result follows. ■ 
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3 GARCH with Gaussian innovations 



The hedging strategies that come out of (2.11 )-( 2.14) are, in general, difficult to compute either explicitly 



or by Monte Carlo methods. Moreover, the interpretation of the values of the resulting local risk- 
minimizing portfolio as an arbitrage free price for h needs of a minimal martingale measure whose 
existence is not always available. 

The approach that we take in this section consists of dropping the physical measure and of carrying 
out the local risk minimization program for a well chosen Girsanov-like equivalent martingale measure; 
we will justify later on that this measure can be used as a legitimate proxy for the physical probability. 
We will implement this program for GARCH models whose innovations are Gaussian, for which no 
minimal martingale measure exists, according to Proposition |2.7| 

The use of a martingale measure for local risk minimization is particularly convenient since the 
formulas that determine the generalized trading strategy are particularly simple and admit a clear 
interpretation. Indeed, it is easy to show that, when written with respect to a martingale measure, the 
local risk-minimizing strategy is determined by 

6 = -^Ek-i [h {ii + auck)] , fc = 1, . . . , T, (3.1) 
Vk^Ekih], fc = 0, . . . , T. (3.2) 

The position on the riskless asset is given by := Vk — ^kSk- Moreover, local risk-minimizing trading 
strategies computed with respect to a martingale measure also minimize |FSc04t Proposition 10.34] 
the so called remaining conditional risk, defined as the process i?^(C°)C) ■= Et[{CT — C**)^], t — 
0, ...,T; this is in general not true outside the martingale setup (see |Sch01[ Proposition 3.1] for a 
counterexample) . 

As we will see in Proposition [33j apart from the computational convenience and the other arguments 
provided above, the chosen equivalent martingale measure has a particular legitimacy since a linear 
Taylor expansion in the drift term of the local risk minimizing value process with respect to this measure 
coincides with the same expansion calculated with respect to the physical measure; consequently, since 
in most cases the drift term is very small, carrying out the risk minimizing program with respect to the 
physical measure or the equivalent martingale one that we introduce below yields virtually the same 
results. 

Consider a GARCH process driven by Gaussian innovations, that is, {ei}ig{i,...,T} ^ IIDN(0, 1). 
Since our intention is carrying out the quadratic hedging program, a challenge at the time of finding 
an equivalent martingale measure consists of making sure that, after the change of measure, we do not 
leave the square-summable category; as we will see in our next theorem this will be ensured by working 
with processes with finite kurtosis. 

Moreover, it is desirable that the innovations do not lose the Gaussian character in the new picture; 
this condition is sometimes imposed as a hypothesis (see, for example, Assumption 2 in jHNOOj ). In 
the next theorem, this is naturally obtained as a consequence of the construction. The proof of the 
following result can be found in the appendix. 

Theorem 3.1 Let (i7, P, J-) be a probability space. Let {sq, si, . . . , st} be a GARCH process determined 
by a recursive re 
J"j :— a{ei, . . . , e 

(i) The process 



by a recursive relation of the type (2.6)-(2.1) and where the innovations {ei}ie{i,...,T} ^ IIDN(0, 1); let 
Ti :— cr(ei, . . . , Ei) be the associated filtration of J- . Then, 



Z„ := n exp ^e, j cxp 1 ^ j , n = 1, . . . , T, 

is a V-martingale. If the drift process {fin}neN is bounded then {Zn}neN is square summable. 
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(ii) Zt defines an equivalent measure Q such that Zt = ^ • 

(iii) The process 

+ n = l,...,r, (3.3) 
forms a IIDN(0, 1) noise with respect to the new probability Q. 

(iv) The log-prices {sq, si, . . . , st} form a martingale with respect to Q and they are fully determined 
by the relations 

Sn = so+o'ie'iH l-CTn^n, (3.4) 

~ '^n (^T^— 1 7 ■ ■ ■ 7 max(p,g) 1 ^n— 1 ; • ■ • ; f^n — g) ■ ('^•^) 



using 



The functions cr„ are the same as cr„ in (2.1) with e„_i, . . . , e„_g written as a function o/e„_i, . . . , e„ 



3.3). If the process {o'nen}ne{i,...,T} is chosen so that it has finite kurtosis with respect to . 



and the drift process {fJ.n}n<£Ti is bounded, then the martingale {sq, si, . . . , st} is square integrable 
with respect to Q. 

(v) The random variables in the process {(Tiei}ig{i t} o-f^ zero mean and uncorrelated with respect to 
Q. 

Remark 3.2 The conclusion in part (iv) on the Q-square integrabihty of the martingale {sq, si, . . . , st} 
is very important since it spells out clearly sufBcient (but not necessary!) conditions, namely, finite 
kurtosis and bounded drift under which we are entitled to use this new measure theoretical representation 
of the problem to compute hedging strategies via local risk minimization. In some market conditions, 
the finiteness of the kurtosis may be a rather restrictive condition but it is well characterized in the 
GARCH context (see |LLMc02[ ILMc02[ ILMc02b| and references therein). 

The condition on the finiteness of the kurtosis can be weakened to requiring the process {cynf^n}n£{i,...,T} 
to belong to L^+'(r2, P, J^), with e > arbitrarily small. This result follows from using in the proof 
(available in the appendix) the fact that the elements of the process {^n}n6{i,....T} do actually belong 



to i''(r2, P, J^), for any q < oo and by replacing the Cauchy-Schwarz inequality in (5.16) by Holder's 
inequality. 

The local risk-minimizing strategy associated to the martingale measure. Given a European 
claim H{St) on the risky asset 5, the martingale measure that we described in the previous theorem 
can be used to come up with a local risk-minimizing strategy by recasting the problem as a hedging 
problem where we consider the log-prices s„ as the risky asset and /i(st) '■= H{g'^v{st)) as the payoff 
function. 

Suppose that the process {<^n^n\n£{Q,...,T} lias finite kurtosis with respect to the physical measure 



P and that the drift is bounded. Part (iv) of Theorem 3.1 guarantees in that situation that the 
log-prices {so,...,st} are square integrable martingales with respect to Q and hence the local risk 
minimization approach to hedging applies in this transformed setup. A straightforward computation 



using the elements in Theorem 3.1 shows that, for any n G {1, . . . , T}, 



£;„_i[s„] = s„_i, £;„_i[(s„ - Sn-iY] = £^„-i[(cr„e„)^] = cr,^^, and var„_i[s„ - s„_il = al 



With these elements, the general local risk-minimizing strategy described in (2.11 )-(2.14) becomes, with 
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the use of this measure: 



Vk ^ Ek[h{sT)], k = 0,...,T, 



£,k — Ek-l 



EfeVfe 



—Ek-i \ekEk[h{sT)] 



= — Ek^i[ekh{sT)], fc = l,...,r, 



(3.6) 
(3.7) 



Lj. = Ct-Co^ Hst) - Vo -J2Usk - Sk-i) = h{sT) - E[h{sT)] -J^^kEk-i [ekHsr)] . (3.6 



fe=i 



fc=i 



The position on the riskless asset is given hy S,'^ := Vk — £,kSk- 

We conclude this section by showing that, since in practice the trend term /i„ is usually very small 
when the time scale is days or weeks^, the value process for the derivative product h obtained by risk 
minimization using the martingale measure that we just introduced and the one computed using the 
physical measure, are very close. We make this explicit in the following statement, whose proof is 
provided in the appendix. 

Proposition 3.3 Let h be a derivative product whose underlying asset is modeled using a GARCH pro- 



cess with finite kurtosis and constant drift. Let Vk be the value process (2.13) of the local risk minim izing 



strategy associated to h computed with the physical probability. Let Vk be the value process (3.6), this 
time computed with respect to the martingale measure introduced in Theorem \3.1\ The linear Taylor 
expansions of Vk and Vk in the drift term coincide. 



3.1 The martingale measure in the price representation and local risk min- 
imization 



The martingale measure presented in Theorem |3.1| is constructed so that, after risk neutralization, log- 
prices become martingales driven by IIDN innovations. Obviously, this feature does not guarantee that 
the prices themselves share this attribute. However, Theorem |3.1| can be easily modified so that an 
equivalent result is available for prices. Indeed, consider a European derivative product with payoff 



function H, whose underlying asset is modeled using a general GARCH process as in (2.6l-(2.7), driven 
by IIDN(0,1) innovations, that is 



log 



Sy-, 



Sn-1 



("'ri — 1 J ■ ■ ■ ^ "'n— max()9,g) ^ ^n— 1 1 ■ • ■ i ^n—q)j 



with {e„}„gp 



where 



and consider 



IIDN(0, 1) and {/i,i}„gN a predictable process. Rewrite (3.9) as 



(3.9) 
(3.10) 

(3.11) 
(3.12) 



an 



n = 1, 



,T. 



■^As an example consider the drift term ^ corresponding to the standard GARCH models with constant drift term 
historically calibrated to the daily log-returns of the following indices between the dates January, 2nd 2007— December 
31st, 2008: Dow Jones Industrial Average: -6.99- lO—*, Nasdaq Composite: -S.SS-IO"", S&P 500: -8.94-10~'', Euronext 
100: -1.1 ■ 10^3. 
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The proof of Theorem |3.1| can be mimicked to show that the measure Q determined by the Radon- 
Nikodym derivative Zt dQ/dP given by 



fe=i 



( 

CTfe 



e/c exp 



2ai 



(3.13) 



is such that the process {e„}„g{i....^T} is an IIDN(0,1) noise with respect to it. Moreover, using 
{^n}ne{i....,T} we have that ( 3.9 )-( 3.10) become 



^ n ; • ■ ■ 1 f^n— max(p,(7) i ^n — 1 7 ■ ■ • ; ^n — q)- 



(3.14) 
(3.15) 



The importance of this modification is that the price process {Sn}ne{o,...,T} forms a martingale with 
respect to Q. Indeed, using the facts that Sn/Sn-i = exp (— ^cr^ + cTnf^n), {en}n6{i,...,T} is an IIDN(0,1) 
noise with respect to Q, and that cr„ is J>i_i-measurable, we obtain that 



En- 



Sn-1 



— En-1 



exp [ --o-„ + cr„e„ 



= exp ( -^(^l ] En-i [exp (cr„e„)] = 1, 



as required. The local risk-minimizing strategy associated to this martingale measure is easy to compute 
and is given by the expressions: 



Vk = Ek[H], k = 0,...,T, 
? _ Sk^i ~ 



H [cxp[ -^cTk +(^k^k ] - 1 



, fc — 1 , . . . , T, 



(3.16) 
(3.17) 



fc=i 



(3.18) 



where E^. := varfe_i(S'j, — S'fc_i) = S'^-i i^'^^'' ^ 1^ The position on the riskless asset is given by 



3.2 



Local risk minimization and the pricing formulas of Duan and Heston- 
Nandi 



The martingale measure Q introduced in (3.13) can be used to provide an alternative risk-minimization 



interpretation to the pricing formula in Duan [D95 , which was introduced using a utility maximization 
argument. The same holds for the formula introduced by Heston and Nandi in [HNOOj : in this case, our 
interpretation is even more valuable for, as far as I know, there is no other mathematical argument that 
supports the use of the formula suggested by the authors, other than the fact that the result can be 
expressed using an extremely convenient closed form expression. Moreover, after the change of measure, 
an additional hypothesis is needed in |HNOO| (see Assumption 2 in page 590) in order to make sure that 
the innovations in the model remain IIDN after risk neutralization. We are able to drop that hypothesis 



since the IIDN character of those innovations is in our case part of the thesis of Theorem 3.1 



As an additional bonus, in both cases the pricing via local risk-minimization comes together with 
an associated hedging strategy that does not exist for the model in | HN00j and, as we mentioned in the 
introduction, has been the object of various conflicting proposals in the case of Duan's model (see |GR98) 
for a discussion). 
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By construction, the self-financing hedging strategy associated to the generahzed local-risk minimiz- 
ing strategy that we propose is also a mean-variance optimal strategy when the corresponding mean 
square hedging error is measured with the martingale measure. Nevertheless, the hedging risk is "felt" 
with the physical probability measure. Since there is no theoretical argument that proves that the op- 
timality of the hedging ratios that we propose with respect to others in the literature survives when we 
change to that measure, we will carry out a numerical study in the next section that seems to indicate 
that this is indeed the case. 



The formula (3.171 that provides the hedging ratios needs in general to be evaluated via a Monte 
Carlo simulation. However, it is worth noticing that it is readily available for any payoff H, which yields 
a competitive advantage with respect to sensitivity methods which usually require the computation of 
derivatives of the option price and may prove to be an extremely convoluted task for exotic derivatives. 
The proof of the following Proposition can be found in the appendix. 

Proposition 3.4 The formulas proposed by Duan f D95^ and Heston and Nandi fHNOO^ for the pricing 
of European options that have underlying risky assets modeled via a NGARCH and an asymmetric 
GARCH process, respectively, coincide with the prices of the local risk-minimizing strategies constructed 



using the martingale measures introduced in (3.13), associated to those two processes. 



3.3 Numerical test of the hedging performance of the local risk minimization 
(LRM) scheme 

We just saw in Proposition |3.4| that the local risk-minimizing approach recovers the pricing formulas 
of Duan [D95] and Heston and Nandi [HNOO and we mentioned that there is a general result (see 
for example |FSc04i Proposition 10.37]) that ensures that the associated self- financing trading strat- 
egy ( 3.16 )-( 3.18) is variance optimal, that is, 



E 



k=l 



< E 



H-Vo-J2^k(.Sk-Sk- 



k=l 



(3.19) 



for any self-financing trading strategy ^i,...,^t with initial value Vq. Given that the hedging risk 
is measure d with the physical probability, a question that needs to be adressed is the validity of the 
inequality (3.19) with respect to that measure, that is, replacing in it the Q-expectation E by the 



P-expectation E. An obvious continuity argument suggests that that statement is going to hold true 
provided that the modified drift term {/i„}„gN is "sufficiently small"; however there is no general result 
that we can invoke and we will hence carry out a numerical experiment to asses the performance we are 
interested in. 

Since there is no proposal for hedging ratios in |HNOO| we will limit our comparison to the hedging 
strategy proposal in Duan's model for the European call option (see |D95t Corollary 2.4]) and the 
standard Black-Scholes model. 

More specifically, we will consider an European call option whose underlying asset has as price 

{Sn}ne{Q,. 

0.2, /3i = 0.7, and A = 0.01. Additionally, the initial value equals Sq — 100. 



T} a realization of the model spelled out in ( 5.21 |-(l5.22 1 with p=q=l, ao = 0.00001, ai 
0.01. Additionally, the initial value equals Sq 



Our goal is comparing the three mean square hedging errors E 



H-Vo 



E 



and E 



BS 



J2k=l ^k^i^k 



Y.k=l^k{Sk 



5fc-l 



where: 



(Vo, {$fc }fcg{i. ... ,T}) i s the self-financing trading strategy associated to the local risk- minimization 
scheme ( 3.16 )-( 3.18) for Duan's model with the parameter values that we just specified. 



GARCH options via local risk minimization 



15 



• {^{F^ ' {Cfe}fce{i,...,T}) is the self-financing trading strategy associated to the Black-Scholes scheme 
for a lognormal model with constant volatility a set equal to the stationary value of the volatility 
of Duan's model under the physical probability, namely a'^ — ao/ {1 ^ ai — f3i). 

• i^Q^ A^k}ke{i,...,T}) is the self-financing trading strategy given by Duan's scheme and specified 
in Corollary 2.4 of |D95) , namely, 



St 

Sk-i 



I-Xt>K 



(3.20) 



where K is the exercise price of the European call option in question. By Proposition 3.4 Vq = V 



The mean square hedging errors are computed by simulating one thousand random price paths, 
using Duan's model with respect to the physical probability (5.21 )-(5.21 ), as the error is computed 
using the P-expectation. Each of these price paths are hedged using the three different schemes listed 
above; the Duan and the local risk-minimization approaches require the estimation via Monte Carlo of 



the hedging ratios (3.20) and (3.17) at each time step. These Q-expectations are computed using each 



time ten thousand different random paths simulated with the risk- neutralized version (5.23)- (5.24) of 
Duan's model. The martingale property is extremely important in this setup consequently, before the 
expectations are computed these paths are modified using the empirical martigale simulation (EMS) 
technique introduced in jDS98llDGS01j which, additionally reduces the variance of the estimation. The 
variance of the estimation of the hedging error is evaluated by repeating randomly one hundred times 
the estimation of the mean square error. 

The results of the numerical test are presented in Table [XT] and in Figure [3T] and show that: 



• Both Duan's and LRM hedging schemes yield a smaller mean square hedging error than the 
Black-Scholes scheme. 

• The Duan and LRM hedging schemes produce mean square hedging errors that are not significantly 
very different, even though Duan performs slightly better than LRM for moneyness bigger than 
one and LRM performs slightly better than Duan for moneyness smaller than one. 

• The variance of the LRM estimation is always smaller than the one obtained using the Duan and 
the Black-Scholes schemes. 



4 Conclusions 

In this paper, we studied the applicability of the pricing/hedging scheme by local risk-minimization 
(LRM) to European options with one underlying asset that is modeled using a GARCH process. The 
main conclusions of the paper are: 

• Local-risk minimizing strategies exist for options with square summable payoffs even though it 
is basically only in the presence of models with bounded innovations that a minimal martingale 
measure is available in order to interpret the subsequent price as an arbitrage-free price. 

• Since the conditions for the existence of a minimal martingale measure are too restrictive, we 
have carried out the LRM scheme with respect to a well-chosen equivalent martingale measure 
that produces the same results as the physical probability up to first order in the Taylor series 
expansion with respect to the drift. 
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Figure 3.1: Mean square hedging errors for European call options using the Black-Scholes scheme (BS), Duan approach, 
and local risk-minimzation (LRM). 
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Hedging Errors. European Call Options. 






















Maturity— 


6 time steps 






Vlattirity— 


.1 time steps 
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G time steps 




Mattuity^ 


21 time steps 




So/A- 




1.04 


1 


0.98 


0.94 


1.04 


1 


1.02 


0.94 


1.04 


1 


0.98 


0.94 


1.04 


1 


0.98 


0.94 


BS 


MSE 
std 


0.0262 
(0.0082) 


0.1681 
(0.0118) 


0.1054 
(0.0122) 


0.0111 
(0.0071) 


0.0991 
(0.0398) 


0.2402 
(0.0237) 


0.1978 
(0.0214) 


0.0457 
(0.0188) 


0.1573 
(0.0440) 


0.3250 
(0.0271) 


0.2924 
(0.0330) 


0.1124 
(0.0327) 


0.2136 
(0.0344) 


0.3803 
(0.0460) 


0.3986 
(0.1380) 


0.1545 
(0.0285) 


Duan 


MSE 
std 


0.0253 
(0.0063) 


0.1558 
(0.0103) 


0.0998 
(0.0106) 


0.0108 
(0.0060) 


0.0894 
(0.0233) 


0.2188 
(0.0243) 


0.1844 
(0.0173) 


0.0419 
(0.0128) 


0.1457 
(0.0290) 


0.2937 
(0.0202) 


0.2675 
(0.0239) 


0.1006 
(0.0204) 


0.1999 
(0.0301) 


0.3465 
(0.0383) 


0.3411 
(0.0428) 


0.1427 
(0.0255) 


LRM 


MSE 
std 


0.0261 
(0.0060) 


0.1544 
(0.0091) 


0.0979 
(0.0092) 


0.0101 
(0.0052) 


0.0915 
(0.0201) 


0.2201 
(0.0221) 


0.1825 
(0.0162) 


0.0409 
(0.0112) 


0.1536 
(0.0240) 


0.2965 
(0.0184) 


0.2673 
(0.0224) 


0.0991 
(0.0183) 


0.2135 
(0.0275) 


0.3545 
(0.0356) 


0.3401 
(0.0374) 


0.1404 
(0.0240) 


std BS 


vs. Duan 


30.1.5% 


14.56% 


15.09% 


18.33 % 


70.81 % 


-2.46 % 


23.69 % 


46.87 % 


51.72 % 


34.15% 


38.07 % 


60.29 % 


14.28 % 


20.10 % 


222.42% 


11.76 % 


std BS 


vs. LRM 


36.66% 


29.67% 


32.60% 


36.53 % 


98.00 % 


7.23 % 


32.09 % 


67.87 % 


83.33 % 


47.28% 


47.32 % 


78.68 % 


25.09 % 


29.21 % 


268.98 % 


18.75 % 


std Duan vs. LRM 


5.00% 


13.18% 


15.21% 


15.38 % 


15.92 % 


9.95 % 


6.79 % 


14.28 % 


20.83 % 


9.78% 


6.69 % 


11.47 % 


9.45 % 


7.58 % 


14.43 % 


6.25 % 



Table 3.1: Mean square hedging errors (MSE) for European call options using the Black-Scholcs scheme (BS), Duan 
approach, and local risk-minimization (LRM). The MSEs are computed by simulating 1,000 random price 
paths, using Duan's model with respect to the physical probability; the Duan and LRM hedging ratios for each 
of those price paths are computed via Monte Carlo using each time 10, 000 random paths simulated with the 
risk-neutralized version of Duan's model to which the EMS modification is subsequently applied. The standard 
deviations (std) are obtained from 100 independent Monte Carlo hedging error estimates. The rows "std A vs. 
B" show the percentage increase in standard deviation experienced in the estimation when using the method 
A instead of B. 



• When this martingale measure is used in the context of the standard models of Duan |D95j 
and Heston and Nandi |HNOO| we recover their pricing formulas, which provides an alternative 
optimal hedging interpretation to the original utility maximization argument that motivated their 
introduction. 

• The hedging strategy associated to the LRM scheme is additionally mean-variance optimal when 
the mean square hedging error is computed using the martingale measure. 

• We have studied numerically if this hedging optimality survives when the error is measured using 
the physical probability by comparing the hedging performances of the standard Black-Scholes 
scheme, Duan's scheme, and local risk minimization. Simulations show that the Duan and LRM 
hedging schemes produce mean square hedging errors that are not significantly very different, 
even though Duan performs slightly better than LRM for moneyness bigger than one and LRM 
performs slightly better than Duan for moneyness smaller than one. The variance of the LRM 
estimation is always smaller than the one obtained using both the Duan and the Black-Scholes 
schemes. 

• The expression that provides the LRM hedging ratios needs in general to be evaluated via a Monte 
Carlo simulation but it is readily available for any payoff H and, unlike other hedging proposals 
based on sensitivity methods, it does not require the computation of derivatives of the option price 
which may prove to be a very convoluted task for exotic derivatives. 



5 Appendix 



5.1 Proof of Proposition 2.1 



The proof of the full statement in Proposition 2.1 is lengthy and convoluted. The reader is encouraged 
to check with [LMcOfbj ILLMc02] , and references therein. In the following lines we will content ourselves 
with checking that condition (2.31 implies the asymptotic weak stationarity of the solutions of the model 



and we will establish (2.4) 
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We start by noting that S„_i[r„] = = -E[r„], var„_i(r„) = E„^i[rl] - £;„_i[r„]2 = a^, and hence 
var(r„) = E [var„_i(r„)] + var (£^„-i [r„]) = E [a^] . 



We now take expectations on both sides of (2.2), that is, 

p q 

i=l 1=1 

taking into account that E[fn] = 0, E[f'f^] = E[a'^-^e'^] = E[a^], and £'[|r„|r„] = 0. We obtain 

E[<jI] =oj+{1+ j^)A{L)E[al] + B{L)E[al], 

where A{L) and B{L) are the polynomials A{z) — 'YTi^i^i'^^ ^ ^i^) — A-^* on the one-step lag 

operator L. Equivalently, 

E[al] = c + [(1 + j')A{L) + BiL)] E[al]. 

This difference equation is stable (see, for instance. Proposition 2.2, page 34 in |H94p . that is, it admits 
an asymptotic solution whenever the roots of the polynomial 



l_(l+^2)A(z)-i?(z)=0, 



(5.1) 



lay outside the unit circle, in which case, expression (2.4) clearly holds. This condition on the roots 
of (5.1) is equivalent to 



(l+7')A(l) + i3(l) < 1, 



(5.2) 



which coincides with ([2^. Indeed, if (1 +72)^(1) + 5(1) > 1, we have that since (1 +7^)^(0) + 5(0) = 
< 1, then (5.1) has necessarily a real root between and 1. Conversely, assume that (5.2) holds and 
that Zq is a root of (5.1) such that \zq\ < 1. Then, 



l = il + -f^)A{zo) + B{zo) 



i=l 



< (1 + 7') ^ a.|zor + E ^^i^oi' ^ (1 + ^')^(i) + ^(1)' 



i=l 



which contradicts our hypothesis. 



5.2 Proof of Theorem O 

(i) We start by proving that 

E[\Zn\] = E[Zn] - 1, for all n = 1, . . . , T. 



(5.3) 



This equality will be needed later on and guarantees that Z„ e L^{il,V,T). Indeed, let p{x) 
i= exp(— x^/2) be the standard normal distribution. Then, 



/ + 00 
dxi • • ■ dxn exp 
-co 



CTl 2cr2 



p{xi) ■ • - exp 



dxi ■ ■ ■ axn-i exp — ^ 



p{xi) ■ ■ - exp - 



A*?!— l2^ri — 1 Mri— 1 
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dx„ exp 



p{Xn) = 1, 



(5.4) 



Given that 

(T„(x„-i, ...,xi) 2a'i{xn-i, ■ ■ ■ ,a;i), 

and that we can repeat this integration procedure n — 1 times more, we conclude that E[Zn\ = 1. We 
now recall that ui, . . . , cr„, /ii, . . . , fin, as well as ei, . . . , e„_i are measurable and hence we can 



write 



En-l[Zn] — En-1 



n *^^p 
.fc=i 



— ( 



e/c exp 



n "^^p 

k=l 



( 

0-fc 



£k exp 



1a4 

2-^ 



Id 

'2 a?. 



exp 



1/4 



exp 



Since (T„ and fin are J>i_i-measurable and is independent from J-n-i, this can be rewritten as (see, 
for example. Proposition A. 2. 5 in |LL08j ) 



En-l[Zn] = Zn-1 CXp 



2 rr2 

^ '^n / J —00 



+00 



exp 



x] dx ^ Zn^l, 



as required. We conclude by showing that Z„ is square integrable for all n — 1, . . . , T whenever the 
drift term is bounded. Indeed, let _B > be such that fin < B, for all n £ N, then 



+00 



^[^n] = / ' ' ' '^^n exp 



-00 
+00 



2A^ia:i _ f4 



p(a;i) • • - exp 



IflnXn fl^^ 



P{Xn) 



dxi ■ ■ ■ dxn~i exp 



2^*1X1 jil 



p{xi) ■ • - exp 



+ 00 



Given that 



dxn exp 



dx^ exp 



(Ti erf 
(Tn{Xn-l, - ■ ■ :Xi) (T^ 

(■^n— 1 ; • ■ • ; -^l ) 



2Mn— l^^ri— 1 Mn — 1 
2 



p(a;„). 



CT„(a;„_i, . . . cr2(a;„_i, ...,Xi) 



p{xn) = exp(^2/^2) < exp(BVw^), (5.5) 



where the inequality follows from the hypothesis (GARCHl) and the bounded character of the drift, 
we can conclude that 



/+00 
dxi--- dXn-l exp 
00 



'i.fiixi ^ fi{ 



p{xi) • • - exp 



CTn-l 



P{Xn- 



Using repeatedly the inequality (5.5) in the previous formula we obtain 

E[Zl\ < exp{nB^/uj^) < +00, 

as required. 



(ii) Zt is by construction non-negative and (5.3) shows that E[Zt] — P{Zt > 0) = 1. This guarantees 
(see, for example. Remarks after Theorem 4.2.1 in |LL08| ) that Q is a probability measure equivalent 
to P. 

(iii) Denote by E the expectations with respect to Q. Then, for any u e M and n E {1, . . . , T}, we will 
prove that 

i^„„i[e™^"] = i?[e™^"] = e-"'/^ (5.6) 
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The first equality in (5.6) together with Proposition A. 2. 2 in |LL08| show that the random variables 
{ei, . . . , er} are independent. The second equality, together with the uniqueness theorem for the charac- 
teristic function of a random variable (see, for instance. Theorem 4.2 in [FF03 ) shows that the random 
variables {e^i, . . . , er} are normally distributed under Q. Indeed, using the Bayes rule for conditional 
expectations and part (i), we have 



£^n-l[-^Te™'^"] — — En-1 



dXn ■ ■ ■ dXT 



exp 



Zn- 

2 

CTn 20-2 



+•00 



dxn exp 

) 

+ C30 

dxn+1 exp 

-co 

Given that all the integrals 



p{x„)e 



V{Xn) ■ ■ ■ exp 



-Er, 



Zn- 



ax 



2a\ 



Zn-l 

p{xt) 



+00 



dxT exp 



2cr2 



p{xt)- 



+00 



dxi exp 



P{xi) = 1, 



the previous expression reduces to 

/'+00 

En-iW'" 



dxn exp 



(5.7) 



(5.8) 



Regarding the second equality in (5.6 1, we compute 
1 



E[e'' 



E[Zt] 

+ CXD 



dxi ■ ■ ■ dxT 



-00 
+00 



exp 



Zq^e \ " "n J 

ai 2af 



dxi exp 



ai 2al 



p{xi) ■ • - exp 

' + 00 



^1TXT ^J'T 



ax 



2a^ 



p{xt) 



Pixi) ■ 



-f 00 



dxn+1 exp - 



dxn exp 

DO 

P{Xn+l) ■ ■ ■ 



-f 00 



dxT exp 



ax 



2a\ 



p{xt)- 



Using again (5.7) and the second equality in (5.8) we easily obtain that 
as required. 



(iv) Expressions (3.4) and (3.5) follow from substituting (3.3) in (2.6) and (2.7). Recall now that 

E{al^^ = E{E,_^{ale% = E{alE{^^\ = E{a1\. 
Hence, by hypothesis (GARCH2), we have that 

E{a1\ = E\a'ief^ < 00. 



(5.9) 
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Using now (5.9), part (i), and Bayes law of conditional probability, we have that 

E[\a,\] = E[a,] = E[Zt<J,] < [E[Z^]] ^ [E[af]] ^ < oo. (5.10) 
Additionally, since by part (iii) the innovations are Gaussian with respect to Q, we have 



E a,M^E 



a.iEi^i [\ei 



= E 



<JiE Wa 



EWa, 



which together with (5.10) implies that i?[|criei|] < oo. This inequality and (3.4) show that s„ € 
L^{n,Q,F). Indeed, 



E[\snW - -B[|so + fTiei 



-CT„e„|] <E[\s^\]+E[\ai€i\ 



E[\antnW < OO. 



Finally, 



En-l[Sn] = En-l[Sn-l + (Jn^n] — Sn-1 + 0-„i?„_i [e„] = S„_i, 



which proves that {soiSi, . . . , st} forms a martingale with respect to Q. Notice that in the last two 
equalities of the previous expression we used the conclusion of point (iii). 

Suppose now that the variables {cri€i}i^{i,...,T} have finite kurtosis with respect to P and that the 
drift process {/injngN is bounded. Then, for each i e {1, . . . ,T} 



E {at4] < oo. 

Then, since E [crfej] = E [crfE^^i [ej]] = 3E [af] , we have that 

E [af] < oo. 



We will proceed by showing first that (5.11) and (5.12) imply that 

E [4] < oo 



(5.11) 



(5.12) 



(5.13) 




i<] = l 



< oo. 



(5.14) 

When the square inside the expectation is expanded, some algebra shows that E [s^] is a finite sum of 
real numbers plus terms that, up to multiplication by finite constants, have the form: 

• E [aiEiUjej] — E [Ei^i [aiSiajej]] = E [aiejajE [e^]] = 0, where we assume, without loss of gener- 
ality, that j < i. 

• E [afef] < oo, by hypothesis (GARCH2). 



• E[afef] <oo, by (5.11). 



• Also by (5.11), the terms of the form 

,2,2„2,21 ^ / TP r„4,41\l/2 / j-n [^4,41 \ 1/2 



Elafe^a]e^]<{Elat4]r'{Elale^]y'<oo. 



(5.15) 
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• E [aiCTjakiJieiejekQ]- This term is also finite because by ( |5.15[ ) 

\E [a,a,a,aie,e,e,ei] \ < {E [aMay^]f' {E [alelafef]^' < oo. 

• Analogous arguments can be used to prove the finiteness of the remaining terms that have the 
form E [afefajc-i], E [afe^ajejakek], E [crfef], E [afefojej], and E [a,ajake,ejek]. 



This argument establishes (5.14). We now use this relation to conclude that {sq, Si, . . . , st} is square 
integrable with respect to Q. Indeed, by part (i) of the theorem, {Zn}ne{i,...T} is a square integrable 
martingale and hence 



E [si] = E [Zrsl] < {E [Z^f' {E [4])^^^ < oo, 

as required. 

(v)Let n e {1, . . . , T}. Then, by part (iii) 

En-l[On^n] = 0'„i?„_i[e„] = (TnE[€n] = 0. 

Now, as £'[cr„?„] = E[En-i[<Tnen]] — 0, the first statement follows. 

Let j € {1, . . . , T} and assume, without loss of generality, that j < n. Then, 



Consequently, 



-E'„_i[cr„e„(Tjej] = CT„(Tjeji?„_i[e„] = (T„o-jeji?[e„] — 0. 
cov((T„e„,crje-,) = £;[cr„e„CTj?j] = E ^;„_i[cr„e„CTje,] = 0. 



(5.16) 



5.3 Proof of Proposition 3.3 



Let /fe(/i) be the function defined by the value process (2.13) with respect to the physical measure, that 
is, 



fk{^J.) :== Ek 



(Tt 



1 - 



A straightforward computation shows that 



fk{0)^Eu[h] and f'M^- J2 ^'^ 

j=k+i 

Consequently, the linear Taylor approximation VJ?*" of Vk is given by 



Vt = Ek[h]-^i Ek 

J=k+l 



(5.17) 



(5.18) 



Let now fk{p) be the value process with respect to the martingale measure Q in Theorem 3.1 Using 
the martingale property of the process Z„ that gives us the Radon-Nikodym derivative dQ/dP we have 



Ek[h] = --^Ek[ZTh] = ^Ek[ZTh]=Ek 



Zj" . 



= Ek 



hexp 



• • • cxp 



CTfc+i 



(5.19) 
(5.20) 



A straightforward computation shows that /fe(0) = /fe(0) and /^(O) = /^(O). Consequently, V^'™ = V^'™, 
as required. I 
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5.4 Proof of Proposition |3.4 

We start with Duan's model which is given by 



log ( = Sn-l = Acr„ - ^cr,^ + (T„e„, 



(5.21) 



ao + J2(^^^l-^^l-^ + T.f^l''n-P {e„}„eN IIDN(0, 1) (5.22) 

1=1 j=i 



where ao > 0, ai,Pj > and J2i=i '^i + S^=i ft < 1 so that second order stationarity is ensured and 
the coefficient A is interpreted as a unit risk premium. Using the notation introduced in (3.9) and (3.12) 



this is a general GARCH model with ^„ = AtT„ — ^cr„ and ^„ ~ Act„. Consequently, by (3.14)-(3.15 ), 
after risk neutralization this process is given by 



log ( -^^j = ^ '"^n-l = -\'^n + (^n^n, {^n}nen ^ nDN(0, 1) 



(5.23) 
(5.24) 



and hence, according to (3.16), its local risk-minimizing price is given by Vb — E[H], which coincides 
with the formula proposed by Duan since the process (5.23)-(5.24| is identical to the one obtained in 



Theorem 2.2 of [D95| out of his locally risk-neutral valuation relationship. 
As to Heston and Nandi, their process is given by 



log 



Sn - Sn-l = A(T„ + (T„e„, {e„}„eN ^ IIDN(0, 1), 

q p 



1=1 3=1 



(5.25) 
(5.26) 



where ao > 0, a^, Pj > and the roots of the polynomial — + caii)^^ * inside the unit 

circle so that second order stationarity is ensured. This time we have that /i„ — Act^, /I„ = (A + |)cr^, 
and 

en :=£«+ (A+-jcr„. (5.27) 



Hence, the risk neutralized version of (5.25)-(5.26l is 



log 



^ , = Sn - s„_i = - l^al + (T„e„, {e„}nGN ^ HDN(0, 1) 



ao + ^ Oil 1 Cn-% 
i=l 



^ + 7i + 2 ) 



2 P 



(5.28) 
(5.29) 



which coincides with the one proposed by Proposition 1 of |HN00j . We emphasize that by Theorem 3.1 



the risk neutralized innovations {e„}„gN are automatically IIDN(0,1) and, unlike in the treatment carried 
out in jHNOO) . no additional assumption is needed. I 
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